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Abstract 

We analyse LEP and PETRA data on single inclusive charged hadron 
cross-sections to establish new sets of Next-to-Leading order Fragmentation 
Functions. Data on hadro-production of large-px hadrons are also used to 
constrain the gluon Fragmentation Function. We carry out a critical compa- 
rison with other NLO parametrizations. 
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1 Introduction 

Inclusive hadron production in various processes gives the possibility to carry 
out quantitative tests of the QCD-improved parton model. The essential property 
of QCD cross- sect ions which allows comparisons between various processes is the 
factorization property. Indeed thanks to the factorization theorem |]I]], the cross- 
sections are written as convolutions of basic building blocks, such as the quark and 
gluon distributions in the incoming hadrons, the hard subprocesses describing the 
large-angle scattering of partons, and the fragmentation functions of quarks and 
gluons into hadrons. The parton distributions and fragmentation functions can be 
defined in a universal way. Therefore these distributions measured in one reaction 
can be used to perform predictions for another reaction. On the other hand, the 
subprocess cross-section is entirely calculable in perturbative QCD, the only free 
parameter being Aqcd- 

Many processes involving fragmentation functions can be related in this manner. 
Among them the best studied are the e"'"e~-annihilation into hadrons e~^e~ -^ hX, 
the observation of hadrons in DIS experiments ip —>■ ihX, and the hadro- and photo- 
production of large-p^ hadrons hahb — ^ hcX and 'yha — > hcX. Clearly quantitative 
studies of these reactions can only be performed if sets of next-to-leading order 
(NLO) fragmentation functions are available. 

While there are several sets of parton distribution functions for the proton and 
the pion, mainly extracted from DIS data, few sets of fragmentation functions exist. 
It is only recently that many precise data on the hadron energy spectrum in e~^e~- 
annihilation became available from LEP experiments, completing those obtained 
at lower energy at DESY, SLAC and TRISTAN. The accuracy of these data is 
remarkable ; they should allow a determination of fragmentation functions with 
a precision which approaches that obtained in the parton distribution functions. 
A first step towards a NLO analysis of fragmentation functions has been done by 
Chiappetta et al. (CGGRW collaboration) [Q], who studied the fragmentation into 
7r°. As LEP data did not exist at that time, these authors also made a detailed 
study of large-p^ 7r° cross-sections, measured by ISR experiments and the UA2 
collaboration, in order to constrain their parametrizations. LEP data was used by 
Binnewies, Kniehl and Kramer (BKK coll.), in association with lower energy results 



of PEP and PETRA, to determine the parton fragmentation into charged pions and 
kaons, and into neutral kaons p, Q. From these resuhs and the measured ratio p/vr, 
they indirectly obtained fragmentation functions for charged particle^. 

We pursue these studies and propose parametrizations^ of the parton fragmen- 
tation functions into charged particles by directly analyzing the corresponding e+e" 
data. The interest of such functions is that they can be used in all reactions in 
which there is no particle identification. On the other hand this set of fragmenta- 
tion functions can be compared with the BKK set, thus offering an estimation of the 
"theoretical error" embedded in the parametrizations and which comes from various 
theoretical assumptions used to perform fits to data. We shall see, in section 3, that 
the gluon fragmentation function is poorly constrained by e~^e~ data. Therefore 
we shall complete our anlysis by that of the hadron-production of large-p^ charged 
hadrons, a reaction which is sensitive to the gluon fragmentation function. 

This paper is organized as follows. In section 2, we present the theoretical 
expressions we use to calculate the single inclusive hadron cross-section in e~^e~- 
annihilation, and we describe the parametrization of the non-perturbative bounda- 
ry conditions associated to the DGLAP evolution equations for the fragmentation 
functions. We also pay attention to an improved choice of the renormalization and 
factorization scales. In section 3, we describe the fits to LEP and PETRA data and 
discuss our results. Then we study, in section 4, the constraints put on the gluon 
fragmentation function by large-p^ hadron production. We conclude in section 5. 

2 Single Inclusive Cross Section in e^e^ 
Annihilation 

In the QCD-improved parton model, the single inclusive cross-section is given 
by the expression 



dz atot dz a^qgJ^ V ^ ' ^tot dv \ y ^ 

(1) 



■''During the completion of this work, two new papers on fragmentation functions apppeared(5| 
We shah briefely comment on them. 

''These fragmentation functions are available on request (guillet@lapp.in2p3.fr). 



The longitudinal variables z = — — - — and v = — ^-5 — measure the fraction of 
the energy ^/Q^, available in the e"'"e~ CMS, carried away by the hadron h and the 
parton a. The fragmentation functions Da{u,M'^) describe the non-perturbative 
transition of parton a into hadron h. Only the dependence on the factorization 
scale M is perturbatively calculable and given by the DGLAP evolution equation 

M'^^^ = T.P^a^D^{M^) (2) 

where indicates a convolution : (/ C?) g){x) = /q dy dv f{y) g{v) 5{yv — x). The 
time-like kernels are expansions in asiM"^) 

pi{M^) = g(M^) pi:^ + (i(^'))' n?^ + ■ ■ • (3) 

In this paper we work at the NLO approximation and when solving (^, we use the 
NLO kernels P^,^ 0' 3- ^'i- (0) ^^^ ^le solved more easily by working in moment- 
space. We use this formalism here as well as a computer code written by P. Nason 

i- 

The hard subprocess cross-section daa/dv describes the production of a parton a 
in e"'"e~-annihilation. (We consider only the sum of the transverse and longitudinal 
cross sections). It is given by an expansion in as(yU^) where /i^ is the renormalization 
scale : 

^(.. ,^ M^ Q=) ^ ^(.. Q^) + ^y)^(v, M^ q^) + . . . (4) 

We use expressions calculated at order as(/i^) ||lO|- For the total cross-section atot 
we also use the 0{as) expression^] : 

^ ^ ^ ^l--(/^^)) • (5) 



Throughout this paper, we work at NLO accuracy and do not take into account 



the NNLO expressions calculated by Rijken and van Neerven [|T2|. Several reasons 
impose this choice. First, the full NNLO correction is not known ; the 3-loop 
DGLAP kernel is missing. Second, we shall use these fragmentation functions in 



^Expressions of cr^°'"" which also contain the Z" contribution may be found in rcf. ||l^, [T^ . 



the calculation of cross section for which only NLO expressions are available (for 
instance large-px inclusive cross sections) ; it appears more coherent to also use 
NLO parametrizations of the fragmentation functions. One must also notice that 
the NLLO corrections to the e~^e~ inclusive cross section are very small when z is not 
too close to zero or one [|12|. In this paper, we study data in the range .12 < 2; < .9, 
and for 2; ~ .9 the experimental errors are larger than the NLLO corrections. 

In order to neglect kinematical higher twists of order m? /p'^ where ?7i is a hadron 
mass and p its momentum, we restrict the 2;-range {z ~ 2|p|/a/Q^) studied by the 
condition z > .12. This condition also allows us to neglect "MLLA" effects [|T3 | 



which show up at small z and are not included in the NLO formalism described in 
this section. 

On the other hand the 0(0^) corrections to the longitudinal cross section da^/dz 



calculated in [|T2| are non-negligible. However we will not use this cross section 
to constrain to gluon fragmentation functions Dgi^ZjQ"^), because data at large z 
{^ ^ -4) is scarce and not very precise. 



Finally let us mention that we work in the MS scheme and that we use massless 
expressions for the DGLAP kernels and the hard cross sections. However we take into 
account the bottom threshold by starting the evolution of the 6-quark fragmentation 
function at M^ = ml. 

The cross-section (|1]) depends on arbitrary scales which show up in the course 
of the theoretical calculations : the renormalization scale /z and the factorization 
scale M. They must be chosen of the order of y/CP, but their precise value cannot 
be determined by any fondamental rule. However several approaches exist which 
propose prescriptions to improve the choice of these scales ; we adopt here the 
"Principle of Minimum Sensitivity" criterion |l^. At the order ^^(/i^) at which 
daa/dv is calculated (LO calculation in asifi"^)), no prescription constrains /i and we 
choose fi = M. Then we fix M using the PMS criterion 



dM^ \dz ) 



= . (6) 

M=Mopt 

We find by a numerical study of (|^) that M^ ^ is much smaller than Q^. In a large 
range in z, we have Mopt — v^Q^/5. However at small values of 2; (^ < .2), Mopt starts 
to increase to values as large as a/Q^. At Mark II fl^ energy {y/Q^ = 29 GeV), 



the PMS criterion does not work for ^ > -3 ; there is no solution of equation (P). 
Because of this result, we shall only study data with \fQ^ > 35 GeV. This leads us 
to discard PEP data. 

As already discussed, the non-perturbative physics of the cross-section (|l]) is 
entirely contained in the fragmentation functions D^iy^ JVP). The perturbative evo- 
lution of these fragmentation functions is given by eq. (0) which must be supple- 
mented by a non-perturbative initial condition. The latter is fixed by the y-behavior 
of the fragmentation functions at an initial scale Q^. As we are interested in the use 
of the fragmentation functions in reactions where the scale M can be of the order of 
a few GeV (hadro- or photoproduction of large-p_L particles), we start the evolution 
at a small value of M, namely M'^ = Ql = 2 GeV^ (M^ = ml for Df,). At this 
value of M , we fix the shape of the fragmentation function by using a simple and 
standard parametrization 



\Pa y"9 



D,{y,Ql)=N,{l-yy 

Du {y, Qt) = {Nu{1 - yf- + Nu{l - yf-) y"- 

Dd+s (y, Qt) = (Nd+sil - yf'^' + iV,+,(l - yf'^^ 

D,(y,Qt)=N,{l-yf^y''^ 

D,(y,mt)=N,{l-y)S''^ . (7) 

Let us say a few words on the theoretical and "experimental" reasons for this 
parametrization. First we notice that the e"'"e~-annihilation cross-section is only 
sensitive to the sum Dd+s = Dd + Dg. So we are not able to determine Dd and Ds 
separately. This degeneracy could be lifted by looking at DIS data, or at data on the 
production of large-p_L particles. However the large -p± production cross-sections, 
at UAl and UA2 energy, are sensitive to the gluon fragmentation function, and very 
little to the (i-quark fragmentation function. Therefore in this paper, we will not 
give separate descriptions of the d- and s-quark fragmentation functions. On the 
other hand, in order to reduce the number of free parameters of the fit, we assume 
that the small- 2; behaviour of the light quarks are the same : a„ = a^ = a^. 

By performing a fit to e^e~ data (see next section), we observed that several 
parameters of the input ([^) are strongly correlated with each other. For instance 
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Ng, ag and Pg are strongly correlated, as well as A^^, N^, Pu, and a^, Nd+s and A^^. 
The 6-paranieters also are strongly correlated with each other. Because of these 
correlations, the fitting procedure is very lengthy and we cannot obtain a positive 
definite error matrix. Therefore we proceed in two steps. First we perform a fit 
with all parameters free. Then, in order to reduce the correlations, we fix some 
parameters to the values obtained in the first fit. This will be discussed in detail in 
the next section. 

We also observed a strong correlation between the gluon parameters and the 
functional shape of the input distribution Dd+s{y,Qo)- A simple form of the type 
A^(l — y)^y°' leads to a gluon fragmentation function which is in disagreement with 
large-p_L charged particle data of UAl |l^. Thus we use more flexible shapes, as 
those of formula (|^), in order to try to decorrelate, as far as possible, gluon and light 
quark parameters. 

In the present paper, we do not intend to determine from data the value of ^■^^■ 



This value might be sensitive to power corrections |11, 17] and such an analysis is 
beyond the scope of this work. Here we use a flxed value of Ai^^ : Ajg^ = 300 
MeV, in agreement with the CTEQ4M [|18l and MRST [|T9l distribution functions. 
One must note that this value is quite in agreement with the new value obtained in 
Ref|^ in a flt to e~^e~ annihilation data. 

3 Analysis of e+e~-annihilation data 

Inclusive charged particle data from CELLO pf at ^/Q^ =35 GeV, from TASSO 
ig at ^/Q^ =AA GeV, from AMY |2^ at ^/Q^ =58 GeV, from DELPHI [|| and 



SLD ||2^ at ^/Q^ =9L2 GeV are used in the flts. To constrain the flavor dependent 
parameters, the uds, c and b flavor enriched samples from ALEPH [|l^, DELPHI 
p5| and OPAL p6| are used. A b enriched sample from DELPHI [^ has also been 



included. 

To take into account the experimental systematics errors, avoiding the treatment 
of correlated errors, the following procedure has been adopted: the normalization 
errors are not included in the x^ evaluation but a normalization for each experi- 
ment is taken as a parameter of the flt allowed to vary within 3 standard deviations 
of the quoted experimental uncertainty. For ALEPH, as we use enriched samples 



in the fit, we do not used the all charged data which has correlated errors [^ . 
The normalization of ALEPH ligth quark sample is kept fixed to 1., while a com- 
mon normalization Nqpal {Ndelphi 94, Ndelphi) is allowed for all OPAL samples 
(DELPHI 94 samples [|3[, DELPHI |25|). 



Using the program Minuit [0, we performed a first fit to data, obtaining a 
reasonable result with x^ = 215 for 217 data points and 25 parameters. In particular 
the data normalization factors are very close to one. As discussed in the preceding 
section, we observe strong correlations between the parameters and the error matrix 
is not positive definite. 

The strong correlation between P^ and Pu, Pd+s and Pd+s leads us to fix the 
diffferences Pu—Pu and Pd+s — Pd+s to the results of the first fit namely Pu — Pu = 1-15 
and Pd+s — Pd+s = 3.6. Similarly, we fix a^, = and a?, = —1.7. The x^ is now 204 
for 217 points and 14 parameters. But the fit is not yet converging and the error 
matrix is not positive definite. The source of this problem lies in the remaining 
strong correlation between A^^ and P^. Fixing also Ni, and Pi, we eventually get a 
convergent fit. Relaxing the condition a^ = also lead to a convergent fit with a 
positive definite error matrix and the following parameters for the non-perturbative 
inputs at Qg = 2 GeV^. (The data normalization factors Aj are also fixed to the 
values obtained in the first fit) 

D,{y,Ql) = 2.11(1 -y)i-3V^i 

Du{y. Ql) = (3.45(1 - y)2-25 - 0.955(1 - yf-') y-^-^' 

Dd+siy. Ql) = (0.95(1 - y)''' + 10.43(1 - y)^-^^^) y-O'^^ 

D,{y,Ql) = 5.15(1 - y)3-8V°-^« 

D,{y,Ql) = 0.897(1 -yf-'y-'-' 

N ALEPH b =1.03 

Naleph c = 1.014 

Nopal = 0.999 ^^ 

Ndelphi9a =1.0006 

Ndelphi = 1.0005 

NsLD = 0.9998 

NcELLO = 1.0055 

Namy = 1.0117 

Ntasso = 0.9835 

The corresponding x^ is 201 for 217 data points and 13 parameters. 

A technical discussion of errors and correlation matrices is given in appendix 1. 



Here we just quote the results of a full analysis (MINOS [^,p8[) on the parameter 
Ng. We obtain a variation of x^ by one unit when increasing (decreasing) the value 
of Ng to 5.25 (1.10), leaving the other parameteres free. These alternative fits are 
studied and used in section 4 where we discuss large-pt data. 

One notices that the gluon input is flat, comparable to the light quark fragmen- 
tation functions. However e"'"e~ inclusive cross sections poorly constrain the gluon 
fragmentation function and the gluon parameters are determined with a large errror. 
The gluon fragmentation function corresponds to an 0{as) correction to the Born 
cross section, and it is quite sensitive to the functional forms of the light quark in- 
puts. Therefore we cannot draw any physical conclusion from this result. However, 
it turns out that the gluon (||) gives a good description of the UAl large-p^ data. 
This point will be developed in section 4. 

It is interesting to compare our fit to the corresponding data by using ratios and 
linear scales which exhibit agreements or disagreement more clearly. The experi- 
mental error is the quadratic sum of statistical and systematic errors. 

Let us start the discussion by the 6-enriched cross sections (Fig. |1],0). The 
agreement is generally good, but one notices at large z some discrepancies between 
the various data sets which might indicate that the systematic errors are underes- 
timated. The light and c-quark theoretical cross sections and data are displayed in 
Fig. §, Fig. I and Fig. §. CELLO and TASSO data (Fig. |) are mainly sensitive 
to the u and c-quark fragmentation functions, whereas light-quark enriched LEP 
sample (ALEPH uds, DELPHI, OPAL uds, DELPHI uds) are sensitive to the d 
and s-quark fragmentation function. Here one also notices some dispersion of the 
LEP light quark data at large z. The all charged data from DELPHI and SLD, also 
included in the fits, are shown in Fig. p. 

In fig. ^ we compare the result of our fit with the BKK fragmentation function 
[0, ^ at Q^ = 10000 GeV^. The two parametrizations are in reasonable agreement 
for z < .6 ; however one must keep in mind that BKK studied pion and kaon cross 
sections, and indirectly all charged cross sections (see below). In fig. ^ we display 
ratios of fragmentation functions on a linear scale in order to compare them more 
closely. 

The M-quark and 6-quark fragmentation functions are harder in BKK. The effect 
could be understood partly as follows: the ratio of proton/pion cross section, inspired 
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from data (fig. 5 of [|29|), has been given the form 0.195 — 1.35(y — 0.35)^ by BKK. 



However, the high y part of this form, not constrained by data, is negative above 
0.7 ; as a resuh their fit to all charged hadrons up to y = 0.8 will have a tendency 
to enlarge some of the quark fragmentation functions in the large-y region. On the 
other hand, the s-quark fragmentation function is constrained by kaon production 
in BKK approach and cannot be directly compared to our parametrization which 
does not distinguish between s-quark and d-quark ; from Fig. ^ one can deduce that 
the ratio {d + s)/{d + s)bkk increases from 1.1 at y ^ .1 to 1.5 at y ~ .8. The 
gluon fragmentation functions have different shapes but, as already discussed, they 
are not very well constrained by e^e~ annihilation data. 

The BKK collaboration also performed a direct fit of ALEPH charged hadron 
data and we can compare the parametrizations given in the thesis of J. Binnewies 
pO| with our results. The ratios of the fragmentation functions are displayed in Fig. 
^ The agreement between the quark distributions is reasonable for -1 < 2; < .6. 
For z > .6, our fragmentation functions are much larger than those of ref. [0. 
One reason for this difference may come from the choice of the data. We choose 
to fit all LEP data. As a result, our cross sections overshoot ALEPH data in the 
large- 2; domain. On the other hand Binnewies' parametrization slightly undershoots 
ALEPH data. Let us also notice that in ref. |5^ ALEPH data is studied only in 
the range z < .8. Therefore we may expect differences with our parametrizations 
for z > .8. 

Finally, one must also keep in mind that we used "optimized" cross sections 
(expression (^) whereas BKK used the scale M^ = Q^. Moreover the bottom 
threshold is set at M^ = Ami by the BKK collaboration whereas we have M^ = ml. 
This may introduce slight differences between the parametrizations. 

The author of Ref. also compared his parametrization with the BKK one and 
observed disagreements for z > .5 . The behaviors of the ratios shown in Fig. ^ 
have similarities with those shown in Ref. [^ for the quark fragmentation functions. 
But one notes a strong disagreement for the gluon fragmentation function which is 
very small for z > .5 in Ref. [^] compared to the BKK parametrization. Such a 
behavior seems to be in contradiction with the large-pt UAl data discussed in the 
next section, where an 83% (54%) contribution to the cross section is due to the 
gluon fragmentation at pt = 5GeV/c {pr = 21.9GeV/c). 
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4 Analysis of large-px production rate 

The determination of the gluon fragmentation function Dg{z,Q'^) from the in- 
clusive e^e~ cross section is not very precise. The gluon contribution to the cross 
section is a NLO effect, and, during the fitting procedure, we found that the gluon 
parameters were quite sensitive to the functional form at Qo chosen for the quark 
distributions. 

A better constraint could come from the longitudinal cross section dai/dz. But 
data at large z, where we want to determine Dgi^ZjQ"^), is rare and not accurate. 
Therefore we turn to hadroproduction of \aige-p± hadron which constrains the frag- 
mentation function in the large-z region [z ~ .7). In hadronic collisions at small 
x± = 2p_L/A/i, the contribution to the cross section involving the gluon fragmenta- 
tion function are large, and the theoretical predictions are quite sensitive to them. 

Data from the UAl collaboration on charge hadrons at large p± (1 GeV < p± < 
20 GeV) precisely cover this kinematical domain. For the theoretical predictions. 



we use the code of ref. [31] which includes NLO corrections. The quark and gluon 



distribution are from MRST-2 |T^ and Aj^ = 300 MeV. The factorization and 
renormalization scales are set equal to p±/2 ; a choice which is dictated by the 
stability of the cross section when the scales vary around p±/2 [^. Our results 
for the region p_L > 5 GeV/c where we trust higher order QCD correction are 
compared with data in fig. [l^. The slope and normalization are well reproduced, 
theory slightly overshooting data. The contribution to the cross section due to the 
gluon fragmentation is important, going from 83% to 54% from pt = 5GeV/c to 
Pt = 21.9GeV/c. 

We can try to improve the agreement between large-p_L data and theory by 
modifying the gluon parameters (^ within the allowed error range. With Ng = 5.25 
on the one hand, and A''^ = 1.1 on the other hand, we obtain two new sets of 
parameters and a x^ value increased by one unity : x^ = 202. 

Dg{y,Ql) = 5.25(1 - y)i-8V°-°=' 

Du{y, Ql) = (3.46(1 - y)2-26 _ 1.4(1 - y)3.31) ^-0.77 

D,Uy, Ql) = (0.959(1 - yr' + 8.86(1 - yy-'') y-^-^' (9) 

D,{y,Ql) = 4.07(1 - yf-^V"-^^ 
D,{y,Ql) = 0.897(1 -i/)3-^V'-^ 
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Da{y,Q'o) 


= i.i{i - yy-''y-'-'' 




Du{y,Ql) 


= (3.32(1 -|/)2-24_ 0.5(1 - 


_^)3.39)^-0.62 


Dd+siy,Ql) 


= (0.93(1 -y)0-94- 0.5(1 - 


- y)'-'') y-'-'' 



(10) 

D,iy,Ql) = 6.37(1 - |/)^-0V°-6' 
D,iy,Ql) = 0.897(1 -y)3-7V^-7 

The corresponding fragmentation functions are compared to those obtained from 

the input (R) in Fig. ITT] and O for Q'^ = 10000 GeV^. We note a change in the 



gluon shape which is important at large z. Unfortunately the large-p_L cross section 
tests the fragmentation functions in the domain 2; ~ .7 where the change is small, 
and we do not expect a substantial modification of the large-p_L cross section. This 
point is verified in Fig. O and O which shows only a very slight improvement of 



the ratio data/theory. Therefore we conclude that the 3 sets of parameters (|^), (|^) 
and ( |To|) are compatible with UAl data. 

5 Conclusion 

In this paper we have performed a NLO analysis of e+e^ annihilation into 
charged hadrons data in order to determine a new set of fragmentation functions. 
Although many LEP data are now available, including samples with enhanced con- 
tent in heavy quarks, we find that it is impossible to obtain a convergent fit when 
using three (or five for the u-quark) free parameters to characterize the input shapes 
of the quark and gluon fragmentation functions. There are strong correlations be- 
tween the parameters and we do not obtain a positive definite error matrix. 

Fixing some of the strongly correlated parameters to values obtained in a fit 
characterized by a good x^? "we succeed in obtaining a positive definite error matrix. 
The overall agreement with data is quite good. Comparing our parametrizations 
of the fragmentation functions with those obtained by the BKK collaboration [^, 
we note several important discrepancies, especially in the large- 2; domain. These 
discrepancies should come from the assumption made by the BKK collaboration on 
the large- 2; behavior of the charged hadron cross sections. A better agreement is 



reached with the results of ref. [3C] in which charged cross sections are also studied. 
But in the large- 2; region, one again notes important discrepancies. One explanation 
of the discrepancy could be the fact that we use sets of data extending to larger 



2;- values than those fitted in |[30[] . 
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This point also emphasizes the necessity to carefully treat the large- 2; region and 
to resum the large logarithms in the theoretical expressions. In this work this is 
done through the optimized scales which should also be used in the calculation of 
other reactions making use of the present parametrization. 

We also test the gluon fragmentation function in large-p^ hadronic collisions. 
We find that 3 sets of fragmentation functions allowed by e^e^ data within one 
standard deviation and differing in the gluon shape lead to very similar predictions 
for large-p_L cross sections at the SppS energy. The overall agreement with UAl 
data is quite good. 

6 Appendix 

Uncertainties on the parameters of the fit @ can be roughly estimated with the 

error matrix using the curvature at the minimum and assuming a parabolic shape 

27| , ||28|| : the results are given on Table |I| where the parabolic errors quoted take 



into account the effects due to parameter correlations. These errors should be taken 
as lower limits on the errors, as some parameters, which may be correlated to the 
ones estimated, have been fixed in the fitting procedure. The results of a full MINOS 
analysis (non-parabolic chisquare) ^7\, 28], following the x^ out of the minimum and 



finding where it corresponds to A^2 = 1 is given on Table ^ 







Table 1: Errors on 


the parameters and correlation 






par. 


error 






correlation with 










"9 


/39 


Nd+s 


Olu 


Pd + S 


Nu 


Nu 


/3„ 


Nd+s 


Nc 


etc 


13c 


^9 


1.49 


0.98 


0.97 


0.11 


-0.66 


0.10 


-0.003 


-0.16 


-0.03 


-0.71 


-0.6 


-0.6 


-0.5 


"9 


0.53 




0.90 


0.04 


-0.75 


0.05 


0.01 


-0.19 


-0.02 


-0.78 


-0.6 


-0.6 


-0.5 


Pa 


0.40 






0.22 


-0.52 


0.18 


-0.02 


-0.12 


-0.05 


-0.58 


-0.5 


-0.5 


-0.5 


Nd+s 


0.16 








0.20 


0.97 


-0.57 


0.57 


-0.49 


0.10 


0.08 


0.12 


-0.03 


Qu 


0.09 










0.16 


-0.35 


0.55 


-0.32 


0.98 


0.29 


0.36 


0.14 


0d+s 


0.09 












-0.62 


0.58 


-0.57 


0.08 


0.09 


0.12 


0.01 
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-0.92 


0.49 
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-0.05 
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0.25 
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0.33 
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Table 2: Uncertainties on the parameters 



parameter 


parabolic error 


negative error 


positive error 


Ng = 2.12 


1.49 


-1.02 


+3.1 


ag = -0.71 


0.53 


-0.51 


+0.70 


i3g = 1.37 


0.40 


-0.38 


+0.51 


N^+s = 0.948 


0.157 


-0.15 


+0.17 


an = -0.684 


0.099 


-0.11 


+0.098 


(3d+s = 0.944 


0.091 


-0.090 


+0.096 


iV„ = 3.435 


1.40 


-1.56 


+1.35 


N^ = -0.956 


2.01 


-1.94 


+2.19 


Pu = 2.26 


0.25 


-0.33 


+0.22 


7Vd+, = 10.43 


1.83 


-1.95 


+ 1.93 


N^ = 5.15 


1.77 


-1.54 


+2.26 


a^ = -0.789 


0.177 


-0.184 


+0.186 


/5c = 3.87 


0.33 


-0.33 


+0.35 
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Figure 1: NLO inclusive charged particle production in e~^e~ —>■ bX collisions at 
^/s = 91.2 GeV with optimized scales and with fragmentation functions obtained 
here (formula |]) compared to data of the ALEPH |l7] and OPAL collaboration 
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Figure 2: NLO inclusive charged particle production in e~^e~ —>■ bX collisions at 
^/s = 91.2 GeV with optimized scales and with fragmentation functions obtained 
here (formula |^) compared to data of the DELPHI collaboration taken in 1994 (top) 
m and in 1991-1993 @. 
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Figure 3: NLO inclusive charged particle production in e~^e -^ u, d, sX collisions 
at ^/s = 91.2 GeV with optimized scales and with fragmentation functions obtained 
here (formula g) compared to data of the ALEPH Q, OPAL O and DELPHI 
l^ collaboration. 
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Figure 4: NLO inclusive charged particle production in e~^e — > cX collisions at 
^/s = 91.2 GeV with optimized scales and with fragmentation functions obtained 
here (formula g) compared to data of the ALEPH Q, OPAL g] and DELPHI 
IM collaboration. 
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Figure 5: NLO inclusive charged particle production in e^e~ —>■ hX collisions with 
optimized scales and with fragmentation functions obtained here (formula §) com- 
pared to data at ^/s = 35 GeV from the CELLO collaboration pO| and at ^/s = 44 



GeV from the TASSO collaboration 121 . 
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Figure 6: NLO inclusive charged particle production in e+e -^ hX collisions at 
^/s = 91.2 GeV with optimized scales and with fragmentation functions obtained 



here (formula H) compared to data of the DELPHI p3| and SLD p4| collaboration. 
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Figure 7: Quark and gluon fragmentation functions obtained in this paper (solid 
line) are compared at the scale Q^ = 10000 GeV^ to those obtained by BKK [Q 
(dotted line). 
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Figure 8: Ratios of the parton fragmentation functions obtained in this paper (for- 
mula D to those obtained by BKK [ffl. The scale is Q^ = iqooO GeV^ 
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Figure 9: Quark fragmentation functions obtained in this paper (solid line) are 
compared at the scale Q^ = 10000 GeV^ to those obtain in [Q using all charged 
(dotted line). 
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Figure 10: NLO inclusive charged particle production in pp collisions at ^/s = 
630 GeV for fi = M = Mp = Pt/'^ with fragmentation functions obtained here 
(formula ^) and structure functions MRS99-2 |T^ compared to data of the UAl 
collaboration [Q. 
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Figure 11: Ratio of the fragmentation functions given by (P) to those of the best fit 
(I) at the scale Q^ = 10000 GeVl 



27 



^ 1.5 

to 

0) 
j3 

Q 
\ 

1? 1,4 

z 

o 



D(Q'= 10000 GeV') low Ng/best fit 



.3 - 



,2 - 



1.1 - 




0.9 



0,1 0,2 0,3 0.4 0,5 0,6 0,7 0,8 0.9 1 

y 



Figure 12: Ratio of the fragmentation functions given by (|T0D to those of the best 
fit d) at the scale Q^ = iqooO GeV^. 
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Figure 13: NLO inclusive charged particle production in pp collisions at ^/s = 630 
GeV for fi = M = Mp = Pt/2 with the set of fragmentation functions obtained 
here (^ by allowing a Ax^ of 1 to the best fit and structure functions MRS99-2 |19| 
compared to data of the UAl collaboration [Q . 



29 



o 



=§1-6 



1.4 



1 UA1 pp^h X/NLLACGG 



M=^=MF=pt/2 



SF MRS99, FF BFGW low N„ 



1.2 - 

1 — 

0.8 - 

0.6 - 

0.4 - 

0.2 - 





10 15 20 25 30 

Pt [GeV/c] 



Figure 14: NLO inclusive charged particle production in pp collisions at -^i = 630 
GeV for n = M = Mp = Pt/2 with the set of fragmentation functions obtained 



here (10) by allowing a Ax of 1 to the best fit and structure functions MRS99-2 



[19| compared to data of the UAl collaboration [^ . 
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